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1.0  INTRODUCTION 


The  Department  of  Defense  and  the  Environmental  Protection  Agency  (EPA)  Environmental 
Monitoring  Systems  Laboratory  are  conducting  a  major  study  of  atmospheric  visibility  in  the 
California  desert.  The  name  of  this  study  is  RESOLVE- Research  on  Operations- Limiting  Visual 
Extinction.  The  supporting  contractor  for  DOD  and  EPA  is  Lockheed-EMSCO,  Las  Vegas,  Nev., 
who  is  responsible  for  much  of  the  field  monitoring  program. 

The  agencies  and  organizations  conducting  RESOLVE  have  been  establishing  an  extensive 
and  sophisticated  monitoring  program  to  provide  the  basic  data  for  the  study.  It  became  apparent, 
however,  that  before  finalizing  the  monitoring  program,  there  was  a  need  to  settle  certain  issues 
through  literature  surveys,  analyses  of  historical  data  sets,  and  error  studies.  Three  specific  issues 
that  needed  to  be  addressed  were  as  follows: 

1.  Determine  if  the  contribution  of  NO2  to  light  extinction  in  the  study  area  is  low 
enough  so  that  NO2  need  not  be  measured  on  a  routine  basis. 

2.  Assess  the  adequacy  of  24-hour  particle  samples  as  compared  to  8-hour  particle  samples 
with  respect  to  characterizing  visibility  during  daylight  hours. 

3.  Investigate  the  propagation  of  potential  RESOLVE  measurement  errors  through  to  the 
final  modeling  results  for  extinction  budgets  and  particle  source  apportionments. 

The  analyses  of  these  three  issues  are  presented  in  Sections  2  through  4,  respectively,  of  this 
report. 


1.1  ORGANIZATION 

The  work  reported  herein  was  performed  under  a  contract  to  Santa  Fe  Research  Corp.  by 
Lockheed-EMSCO.  Santa  Fe  Research  Corp.  had  the  primary  research  responsibility  for  the 
analysis  of  NO2  extinction  contributions  (Section  2)  and  the  study  of  8-hour  versus  24-hour 
particle  sampling  (Section  3).  The  error  propagation  analysis  (Section  4)  was  conducted  via  a 
subcontract  to  Desert  Research  Institute.  Two  consultants,  Alan  Waggoner  and  Raymond  Weiss, 
provided  the  data  for  Section  3  as  well  as  general  assistance  on  the  other  aspects  of  the  project. 


1.2  SUMMARY  OF  CONCLUSIONS  AND  RECOMMENDATIONS 

The  issue  of  the  contribution  of  NO2  to  extinction  is  investigated  in  this  report  using  three 
data  bases:  measurements  made  from  aircraft,  taken  as  part  of  the  California  Air  Resources  Board 
(CARB)  Los  Angeles/desert  transport  study,  routine  historical  data  for  NO2  and  visual  range,  and 
a  small  amount  of  data  acquired  with  high-resolution  research-type  NO2  monitoring.  All  three  data 
bases  suggest  that,  on  the  average,  fractional  NO2  contributions  to  total  extinction  in  the 
California  desert  are  generally  less  than  5%,  maybe  substantially  less  than  5%.  Accordingly,  it 
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does  not  appear  necessary  to  monitor  NO2  on  a  routine  basis  as  part  of  RESOLVE.  Nevertheless, 
to  become  more  certain  regarding  the  role  of  NO2,  two  further  steps  are  recommended.  First, 
when  hourly  NO2  measurements  are  compiled  as  part  of  the  data  base  for  the  California  Desert 
Air  Working  Group  (CDAWG),  the  measurements  should  be  examined  to  study  NO2  contributions 
under  worst-case  conditions.  Second,  the  RESOLVE  “special  intensive  studies”  should  include 
some  NO2  monitoring  with  high-sensitivity  instruments  at  the  Naval  Weapons  Center  (NWC),  China 
Lake,  and  Edwards  AFB,  Calif. 

A  study  of  historical  nephelometry  data  for  the  California  desert  reveals  that  significant 
diurnal  variations  exist  in  fine  particle  concentrations.  These  variations  do  not  necessarily  imply, 
however,  that  source  mixtures  for  fine  particles  tend  to  vary  significantly  over  the  day.  In  fact, 
a  high  correlation  between  8-hour  (0800  to  1600  PST)  and  24-hour  nephelometer  averages 
indicates  that  daytime  and  24-hour  particulate  concentrations  are  not  uncoupled.  Furthermore, 
an  examination  of  available  particulate  composition  data  with  diurnal  resolution  suggests  that 
daytime  and  nighttime  particulate  source  mixtures  are  the  same  in  an  average  sense.  Because  of 
these  findings,  because  of  procedural  difficulties  and  extra  expenses  associated  with  short-term 
particle  sampling,  and  because  of  some  interest  in  nighttime  as  well  as  daytime  visibility  for  DOD 
operations,  it  seems  prudent  to  begin  the  routine  RESOLVE  program  with  24-hour  particle 
sampling.  It  is  recommended,  however,  that  further  study  of  this  issue  be  conducted.  The  drum 
sampler  data  from  RESOLVE  should  be  analyzed  regarding  daytime/nighttime  differences  in 
particle  contributions  from  sulfates  and  soil.  Moreover,  as  part  of  the  special  intensive  studies  for 
RESOLVE,  the  particle  samplers  should  be  run  on  a  daytime  and  nighttime  basis. 

Least  squares  solutions  to  the  mass  balance  and  scattering  budget  receptor  models  that 
provide  uncertainty  estimates  for  model  results  are  derived  herein.  Applications  of  these  models 
to  randomized  data  derived  from  known  source  contributions  and  known  scattering  efficiencies 
show  that  the  calculated  uncertainties  are  good  estimates  of  the  real  uncertainties,  even  when  the 
measurement  uncertainty  levels  are  as  high  as  30%.  To  characterize  and  minimize  the  uncertainty 
of  the  interpretive  analysis  results,  the  RESOLVE  program  should  document  and  minimize  the 
measurement  uncertainties  and  maximize  the  number  of  sampling  cases. 


2.0  CONTRIBUTION  OF  NO2  TO  LIGHT  EXTINCTION 
LEVELS  IN  THE  CALIFORNIA  DESERT 

Total  light  extinction  consists  of  a  simple  sum  of  four  terms:  light  scattering  by  particles 
(usually  dominated  by  sulfates  and  other  fine  aerosols),  light  absorption  by  particles  (dominated 
by  elemental  carbon),  light  scattering  by  air  molecules  (natural  Rayleigh  scatter),  and  light 
absorption  by  gases  (basically  all  from  NO2).  A  preliminary  working  hypothesis  of  the  RESOLVE 
study  is  that  light  absorption  by  NO2  is  an  insignificant  component  of  regional  extinction  in  the 
California  desert.  This  hypothesis  is  based  on  the  following  argument  In  the  Los  Angeles  Basin 
and  San  Joaquin  Valley,  the  major  source  areas  upwind  of  the  California  desert,  the  average 
contribution  of  NO2  to  total  extinction  is  about  8%.*  We  would  expect,  however,  relatively 
lesser  NO2  contribution  in  the  desert.  One  basic  reason  is  that  transported  particles  have  slow 

*  A  table  of  published  results  regarding  NOj  extinction  contributions  ir  the  source  areas  is  as  follows: 

Los  Angeles  Basin  -  -  -  8%  2  sites  (Reference  1) 

10%  1  site  (Reference  2) 

5-13%  9  sites  (Reference  3) 

San  Joaquin  Valley  -  -  -  7-8%  3  sites  (Reference  3) 


_*w.  _ «-* 
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removal  rates,  while  transported  NO2  is  continuously  depleted  by  oxidation  to  nitric  acid, 
conversion  to  particles,  and  deposition  to  surfaces.  Furthermore,  as  the  particle  and  NO2  pollution 
is  diluted  during  transport,  the  background  Rayleigh  scatter  term  becomes  relatively  more  signifi¬ 
cant,  reducing  the  overall  percentage  contribution  from  both  particles  and  NO2. 

Because  of  the  above  hypothesis,  the  acquisition  of  NO2  data  on  a  routine  basis  is  not 
currently  planned  as  part  of  RESOLVE.  It  is  very  important,  nevertheless,  to  check  the 
assumption  regarding  NO2  levels  in  the  study  area.  The  purpose  of  this  chapter  is  to  estimate 
NO2  contributions  to  extinction  in  the  California  desert  by  analyzing  currently  available  data. 


2.1  CARB  LOS  ANGELES/DESERT  TRANSPORT  STUDY 


No  data  base  exists  that  can  provide  a  definitive  answer  regarding  NO2  extinction  contribu¬ 
tions  in  non-urban  parts  of  the  California  desert.  The  best  data  available  come  from  an  aircraft 
survey  by  Meteorology  Research,  Inc.  (MRI),  Altadena,  Calif.,  conducted  as  part  of  the  Los 
Angeles/desert  transport  study  sponsored  by  CARB  (Reference  4).  The  MRI  survey  may,  however, 
overestimate  typical  NO2  contributions  in  the  desert  because  the  survey  was  designed  to  study 
transport  phenomena  from  Los  Angeles  rather  than  remote  desert  conditions. 

Table  2-1  summarizes  average  particle  scattering  (bSp)  and  average  NO2  concentrations  for 
desert  flight  paths  during  the  survey.  The  percent  contribution  of  NO2  extinction  (bjqQj 
or  bag)  to  total  extinction,  bag/(bSg  +  bSp  +  bag  +  bap),  is  calculated  from  these  data  as 
follows: 

bag  =  NO2  absorption  at  550  nm  (nanometers)  =  0.033- [NO2],  with  units  of  bag  in 
10~*m_l  and  units  of  [NO2]  in  pphm  (parts  per  hundred  million)  as  listed  in 
Table  2-1. 

bSg  =  Rayleigh  scatter  =  0.12  lO^nr*. 

bSp  =  particle  scattering  =  as  listed  in  Table  2-1. 

bap  =  particle  absorption  =  assumed  to  be  10%  of  bSp. 


TABLE  2-1.  Aircraft  NO2  and  BSCAT  Data  by  MRI, 
July  —  August  1981. 


Date 

Time 

Site 

bw  dO'^m'1) 

NO2,  pphm 

NOj  contribution  to  total 
estimated  extinction, 

7/09 

1622 

Cajon  Pass 

0.46 

1.0 

5 

7/09 

1655 

Cajon  Pass 

0.41 

0.8 

4 

7/09 

1720 

Cajon  Pass 

066 

1.2 

4 

7/14 

1750 

Palmdale 

0.57 

1.0 

4 

7/14 

1 908 

Victorville 

0.87 

1.8 

5 

7/15 

0913 

Palmdale 

0  29 

0.3 

2 

7/15 

1008 

Bars  tow 

0.38 

0.7 

4 

7/18 

1800 

Palmdale 

1.13 

1.5 

4 

7/27 

1645 

Victorville 

0.59 

09 

4 

7/30 

1530 

Victorville 

0.33 

0.3 

2 

a  NOo  fractional  contribution  *  - — 

bifl 


t^p  •  bjg  *  b,p,  where  t  *  scattering. 


a  -  absorption,  g  -  gas,  and  p  -  particles. 


\  • 

rj. 
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Table  2-1  shows  that  the  contribution  of  NOj  to  total  extinction  during  the  MRI  surveys 
averaged  2  to  5%.  As  noted  above,  these  contributions  may  be  higher  than  those  typically  found 
in  remote  parts  of  the  desert  not  under  the  direct  influence  of  transport  from  Los  Angeles. 

2.2  ROUTINE  LONG-TERM  DATA  BASES 

Another  way  of  estimating  the  fraction  of  extinction  caused  by  NO2  in  the  California  desert 
is  to  analyze  the  historical  data  bases  for  routine  NO2  monitoring  and  airport  visibility.  Trijonis 
(Reference  S)  has  published  an  isopleth  map  of  median  airport  visibility  levels  throughout 
California,  including  the  desert.  Median  (1300)  visual  ranges,  taken  from  that  map,  can  he 
converted  into  estimates  of  total  extinction  using  a  Koschmeider  constant  of  3.0  (b  =  3.0/V),  the 
appropriate  constant  for  airport  visibility  data  (Reference  3). 

Long-term  continuous  NCh  data  are  available  from  CARB  and,  to  a  lesser  extent,  from 
Southern  California  Edison  (SCE)  Co.  Values  of  NO2  at  1300  were  compiled  from  the  CARB 
data  so  that  they  would  correspond  to  the  1300  airport  visibility  data  (only  annual  means  of  all 
hours  were  readily  available  for  the  SCE  NO2  data).  The  CARB  NO2  data  are  rather  coarse  for 
desert  conditions  in  the  sense  that  values  are  reported  only  to  the  nearest  pphm.  In  determining 
median  NO2  from  the  data,  cumulative  frequency  distributions  were  calculated  assuming  0  pphm 
recordings  meant  concentrations  below  0.5  pphm,  I  pphm  recordings  meant  concentrations  between 
0.5  pphm  and  1.5  pphm,  etc.  These  frequency  distribitions  appeared  generally  well-behaved  when 
plotted  on  log-normal  paper;  the  medians  were  estimated  by  interpolations  of  lines  fit  to  the 
log-normal  plots.  The  median  NO2  concentrations  were  translated  into  NO2  extinction  levels  using 
the  formula  discussed  in  the  ;  revious  section. 

The  results  of  the  analysis  are  presented  in  Table  2-2.  These  results,  although  more  variable 
than  those  in  Table  2-1,  are  consistent  with  Table  2-1  in  suggesting  that  the  fraction  of 
extinction  caused  by  NO2  in  the  desert  averages  about  4%.  It  should  be  noted  that  we  strongly 
suspect  that  these  findings  represent  overestimates  of  actual  NO2  contribution  to  extinction  in 
non-urban  areas  of  the  desert.  The  reason  is  that  most  of  the  sites  in  Table  2-2  involve  a 
comparison  of  NO2  in  towns  or  small  urbanized  areas  with  a  long  path  measurement  of  visual 
range  extending  over  rural  areas.  The  NO2  concentrations  in  the  small  urbanized  areas  are 

probably  greater  than  the  average  NO2  concentrations  over  the  site  path.  Finally,  it  should  be 

noted  that  the  variations  from  site  to  site  in  Table  2-2  do  not  seem  readily  explainable  and  may 
be  due  to  the  inadequacies  of  this  somewhat  crude  analysis. 

We  had  also  hoped  to  analyze  these  data  sets  for  fractional  NO2  contributions  to  extinction 
under  individual  “worst-case”  conditions.  This  analysis  would  require  an  airport  located  near  a 
non-urban  desert  monitoring  site.  It  turned  out  that  no  such  pairing  of  sites  existed.  In  the 
future,  we  hope  to  obtain  bSp  and  NO2  data  from  one  or  two  remote  SCE  sites  that  can  be 

studied  regarding  worst-case  episodes.  Such  data  will  become  available  when  a  data  base  compila¬ 

tion  study  now  being  conducted  by  CDAWG  is  completed. 

2.3  HIGH-SENSITIVITY  NO2  MONITORS 

The  above  analyses  are  based  on  data  from  “off-the-shelf’  NO2  monitors.  Because  NCb  con¬ 
centrations  are  very  low  in  the  desert,  the  observations  generally  occur  near  the  detection  thresh- 
holds  of  the  instruments.  We  have  found  one  NO2  data  set  for  the  Mojave  Desert  taken  with  a 
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high-sensitivity,  low-zero  uncertainty  NO2  monitor  (Reference  6).  For  these  measurements,  taken 
during  September  1983,  NO2  concentrations  were  always  below  1  pphm,  and  generally  below  0.5 
pphm.  Assuming  a  typical  total  extinction  level  of  0.4  lCT^m-',  NO2  concentrations  of  0.5  pphm 
would  represent  4%  of  total  extinction. 


TABLE  2-2.  Median  Midday  NO2  and  Visibility  Data  for  Desert  Sites. 


City 

(population) 

Median  (1300)  NOj, 
pphm 

Median  (1300) 
annual 

visual  range,  mi. 

NOj  contribution 
to  extinction,  % 

CARB  DATA 

Victorville 

(11,000) 

0.6* 

1 

33* 

4f 

Barstow 

(17.000) 

C.7 

36 

5 

Lancaster 

(35,000) 

0.5 

30 

3 

Indio 

(15.000) 

0.7 

30 

4 

Palm  Springs 
(21,000) 

1.0 

24 

4 

Trona 

(1.500) 

0.6 

65 

7 

Palo  Varda 
(07) 

0.7 

60 

8 

SCE  DATA 

Coolwat  ar 
(2,0007) 

1.1rf 

36 

7 

Lucarna  Valley 
(0?) 

0.3 

35 

54 

*  Provided  by  CARB. 

*  (Reference  5) 

c  Uiing  a  Koschmeider  constant  of  3.0  (5%  contrast  level)  to  convert  airport 
median  visibility  into  extinctions,  end  using  the  NOj  extinction  coefficient  given 
in  the  text. 

d  Provided  by  SC.E, 


2.4  CONCLUSIONS 

All  of  the  analyses  conducted  in  this  chapter  suggest  that  fractional  NO2  contributions  to 
total  extinction  in  the  California  desert  are  generally  less  than  5%,  maybe  significantly  less  than 
5%,  given  that  the  techniques  tended  to  yield  overestimates.  Accordingly,  it  does  not  appear 
necessary  to  monitor  NO2  on  a  routine  basis  as  part  of  RESOLVE.  Nevertheless,  to  become  more 
certain  regarding  the  role  of  NO2,  two  further  steps  are  recommended.  First,  when  the  CDAWG 
data  base  becomes  available,  some  of  the  SCE  data  for  NO2  and  bSp  should  be  examined  to 
study  NO2  contributions  under  episode  conditions.  Second,  the  RESOLVE  special  intensive 
studies  should  include  some  NO2  monitoring  with  high-sensitivity  instruments  at  NWC  and 
Edwards  AFB. 
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3.0  ADEQUACY  OF  24-HOUR  PARTICLE  SAMPLING  SCHEDULE 


With  respect  to  the  current  RESOLVE  monitoring  plan,  there  is  a  question  as  to  whether 
the  24-hour  fine  particle  samples  are  adequate  to  address  daytime  visibility  problems.  Specifically, 
concern  exists  whether  daytime  particle  source  mixtures  might  differ  significantly  from  24-hour 
average  source  mixtures.  The  purpose  of  this  section  is  to  evaluate  the  adequacy  of  24-hour 
particle  samples  as  compared  to  8-hour  daytime  particle  samples.  This  issue  will  be  addressed  by 
conducting  three  types  of  analyses:  assessing  the  importance  of  and  causes  for  diurnal  variations 
in  existing  nephelometry  data  (Section  3.1),  examining  the  predictability  of  8-hour  average  bsp 
(fine  particle  scattering)  from  24-hour  average  bsp  (Section  3.2),  and  studying  daytime/nighttime 
particulate  composition  data  published  in  Ouimette’s  Ph.D.  thesis  (Section  3.3)  (Reference  7). 

Each  of  the  three  planned  analyses  has  significant  limitations.  The  first  two  are  based  on  a 
large  volume  (several  years)  of  nephelometry  data,  but  they  do  not  explicitly  address  the  problem 
of  diurnal  variations  in  particulate  chemical  composition  (i.e.,  particulate  source  mixture).  The  third 
analysis  examines  actual  chemical  composition  data,  but  it  is  based  on  only  a  few  days  of 
sampling  in  one  season.  Combining  all  three  analyses,  however,  should  provide  fairly  sound  recom¬ 
mendations  for  the  RESOLVE  sampling  program  (Section  3.4). 


3.1  DIURNAL  VARIATIONS  IN  NEPHELOMETRY  DATA 

One  way  of  shedding  light  on  the  issue  of  day/night  variations  in  fine  particle  source 
mixtures  is  to  examine  the  diurnal  patterns  in  the  historical  nephelometry  (bsp,  particle  scattering) 
data  base  for  the  study  area.  Are  the  diurnal  patterns  in  bsp  (i.e.,  fine  p article  concentrations) 
important  or  insignificant?  Are  the  diurnal  patterns  likely  associated  with  significant  shifts  in 
particle  sources? 

Figures  3-1  through  3-5  plot  annual  average  diurnal  patterns  for  five  DOD  nephelometer 
sites  in  the  California  desert.  In  order  to  avoid  seasonal  biases  in  the  annual  plots,  the  analysis 
is  limited  only  to  years  for  which  complete  data  were  available  at  the  time  of  this  study.  The 
analysis  is  further  restricted  to  days  providing  at  least  18  hours  of  valid  data  (nearly  all  days 
meet  this  restriction).  Diurnal  patterns  for  each  site  on  a  monthly  basis  are  presented  in  the 
appendix  to  this  report. 

Each  appendix  figure  contains  two  lines,  representing  two  ways  of  calculating  the  diurnal 
patterns.  The  solid  line  simply  represents  the  annual  average  value  for  each  hour;  it  is  measured 
by  the  ordinate  axis  labeled  “bsp.”  The  dashed  line  represents  an  average,  over  the  year,  of  the 
fractional  change  away  from  the  mean  value  each  day.  The  latter  index,  measured  on  the  left 
ordinate  axis  labeled  “(bsp  -  bave)/bave,”  is  an  indication  that  weights  days  of  low  concentra¬ 
tions  on  an  equal  basis  with  days  of  high  concentrations  in  calculating  the  diurnal  patterns. 

With  the  exception  of  the  Death  Valley  results,*  Figures  3-1  through  3-5  reveal  strong 
diurnal  patterns  in  fine  particle  concentrations  (i.e.,  bsp),  with  a  distinct  minimum  in  the  early 


*  There  are  several  alternate  hypothecs  regarding  the  lack  of  a  significant  diurnal  pattern  at  Death  Valley. 
One  conjecture  is  that  transported  pollution  may  be  homogenized  on  a  very  large  scale  by  the  time  it  reaches 
the  distant  site  at  Death  Valley,  so  that  diurnal  dispersion  variations  do  not  necessarily  produce  diurnal  changes 
in  concentrations. 


afternoon.  One  possible  explanation  is  that  local  source  emissions  (e.g.,  from  woodburning  or 
spaceheating  in  the  winter)  might  accumulate  during  the  evening  to  produce  the  observed  diurnal 
pattern.  This  might  imply  a  distinct  day/night  difference  in  particle  source  composition.  However, 
another  possible  explanation  is  that  fine  particle  concentration  elevation  in  the  late  afternoon 
and  evening  is  due  to  transport  from  upwind  air  basins  and  that  the  daytime  minimum  simply 
reflects  the  dilution  of  this  air  mass  with  increased  mixing  height  during  the  afternoon.  This  latter 
explanation,  which  seems  supported  by  the  fact  that  the  diurnal  pattern  is  greatest  during  the 
si miner  (see  Appendix),  would  be  consistent  with  an  equivalent  source  mixture  during  day  and 
night.  In  any  case,  the  strong  diurnal  patterns  do  suggest  that  the  issue  of  24-hour  versus  8-hour 
sampling  must  be  given  careful  scrutiny. 


O  6  12  18  24 

Hour 


FIGURE  3-1.  Diurnal  Pattern  of  Particle  Scattering  at  China  Lake  Range  Trailer. 
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3.2  CORRELATION  OF  8-HOUR  VERSUS  24-HOUR  NEPHELOMETRY  DATA 

For  further  analysis,  the  nephelometry  data  at  the  five  California  desert  sites  were  organized 
into  three  8-hour  averages  as  well  as  a  single  24-hour  average  each  day.  The  main  purpose  was  to 
correlate  the  24-hour  average  with  the  daytime  (0800-1600)  8-hour  average  (see  lower  righthand 
comers  of  Figures  3-6  through  3-10).  In  addition,  we  used  these  data  sets  to  plot  the  diurnal 
patterns  in  terms  of  8-hour  averages  (upper  lefthand  comers  of  Figures  3-6  through  3-10),  the 
24-hour  average  seasonal  patterns  by  quarters  of  the  year  (upper  righthand  comers  of  the  figures), 
and  the  8-hour  diurnal  patterns  by  quarters  (lower  lefthand  comers  of  the  figures). 


DIURNAL  PATTERN 


SEASONAL  PATTERN, 
24-H0UR  AVERAGE 


DIURNAL  PATTERN 
\  (BY  SEASON) 


CORRELATION  COEFFICIENT: 

BSCAT  24- HR  vs.  BSCAT  8  HR 

ANNUAL  .93 

QUARTER  1  .95 
QUARTER  2  .90 
QUARTER  3  .92 
QUARTER  4  .94 


FIGURE  3-6.  Diurnal  Pat 
far  6-Hour  and  2444 our 


at  China  Lake  R 
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FIGURE  3-8.  Diurnal  Panama.  Saaaonal  Variation*,  and  Comlationt 
lor  8+1  our  and  24-Hour  Partial*  Scattaring  at  Laural  Mountain 


800  1600  2400 


CORRELATION  COEF,  ICIENT: 

8SCAT  24-HR  vs.  BSCAT  8-NR 

ANNUAL  .87 

QUARTER  1  .90 
QUARTER  2  .91 
QUARTER  3  .82 
QUARTER  4  .91 


FIGURE  3-9.  Diurnal  Patterns,  Saaaonal  Variation,  and  Contention 
for  8-Hour  and  24-Hour  Partida  Scattering  at  Randsburg  teaalv 


BSCAT  (10‘V)  BSCAT  <10’ 
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CORRELATION  COEFFICIENT: 

BSCAT  24-HR  vs.  BSCAT  8 -MR 

ANNUAL  .96 

QUARTER  1  .97 
OUARTER  2  .95 
"CARTER  3  .94 
CHARTER  4  .95 


FIGURE  3-10.  Diurnal  Pattern*,  Saannal  Variation*,  and  Corralation* 
for  0-Hour  and  24-How  Parttda  Scattering  at  Oaath  VaNay. 
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Table  3-1  summarizes  the  correlations  between  the  daily  24-hour  average  and  the  daytime 
8-hour  averages  for  each  of  the  five  sites.  The  correlations  are  presented  for  the  entire  data  set 
(annual)  as  well  as  for  the  data  disaggregated  by  seasonal  quarter.  Table  3-1  shows  that  the 
correlations  are  uniformly  high,  generally  between  about  0.90  and  0.95.  This  means  that  24-hour 
average  b$p  can  serve  as  a  fairly  precise  predictor  of  8-hour  daytime  bSp.  It  further  implies  that 
the  daytime  and  24-hour  variations  in  One  particle  concentrations  are  certainly  nor  uncoupled.  This 
is  the  first  piece  of  evidence  hinting  that  24-hour  particle  sampling  may  be  adequate  for  character¬ 
izing  the  8-hour  daytime  period.  Specifically,  even  though  the  8-hour  daytime  averages  are 
generally  lower  than  the  24-hour  averages  (as  demonstrated  by  the  diurnal  patterns),  the  8-hour 
averages  appear  to  be  consistently  related  to  the  24-hour  averages  (as  indicated  by  the  high 
correlations).  This  coupling  further  suggests,  but  does  not  prove,  similar  particle  source  compo¬ 
sitions  for  the  8-hour  and  24-hour  averages. 


TABLE  3-1.  Correlation  Coefficients  Between  Daily  24-Hour 
end  Daytime  8-Hour  Nephelometry  Averages. 


Area 

Annual 

Quarterly 

1 

2 

3 

4 

China  Lake  Range  Trailer 

0.93 

0.95 

0.90 

0.92 

0.94 

B  Mountain 

0.94 

0.95 

0.89 

0.93 

0.94 

Laurel  Mountain 

0.90 

0.85 

0.85 

0.86 

0.94 

Randfburg  Wash 

0.87 

0.90 

0.91 

0.82 

0.91 

Death  Valley 

0.96 

0.97 

0.95 

0.94 

0.95 

3.3  DAY/NIGHT  PARTICULATE  COMPOSITION  DATA 

Currently,  there  appears  to  be  only  one  data  base  available  for  the  California  desert  with 
day/night  resolution  of  aerosol  chemical  composition.  This  is  the  fine  particle  data  set  published 
by  Ouimette  (Reference  7).  Table  3-2  lists  Ouimette’s  monitoring  periods  that  involved  daytime 
samples  immediately  prior  to  or  followed  by  nighttime  samples.  The  table  also  indicates  a  plotting 
notation  to  distinguish  points  in  subsequent  figures. 

Ouimette’s  data  base  contains  seven  daytime  samples  juxtaposed  with  six  nighttime  samples. 
We  have  calculated  an  overall  average  mass  balance  for  these  samples  as  shown  in  Table  3-3.  In 
agreement  with  the  previously  discussed  nephelometry  data,  Table  3-3  shows  that  the  total  fine 
particle  mass  concentrations  are  lower  during  the  day  than  during  the  night.  A  mysterious  and 
as  yet  unexplainable  aspect  of  the  data,  however,  is  the  rather  large  “unaccounted  for”  fine 
particle  mass,  which  is  nearly  identical  in  absolute  magnitude  (5.3  or  5.5  pgjm 3)  for  the  day 
and  night  averages.  Expressing  the  mass  balance  in  terms  of  “percent  of  accounted  for  fine  mass” 
(see  values  in  parentheses)  indicates  that  for  Ouimette’s  overall  average  data,  the  particulate 
chemical  composition  is  essentially  identical  during  day  and  night.  This  fact  is  direct  evidence 
that  the  particle  source  mixture  does  not  change  radically  from  day  to  night,  at  least  on  the 
average. 
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TABLE  3-2.  Sampling  Ptriods  for  Ouimatte't  Day/Night  Particulate 
Competition  Data. 


Typical  templing  conditions 

Midnight  Noon  Midnight  Noon 

* - d7.ZZTi2j.lnZ2 - ^ 


Dayi  and  timat  used: 


Type 

Day 

Start  time 

Sample  length,  hrs. 

O  Day  sample  followed 

8-06 

0900 

6 

by  night  sample 

1900 

11 

8-07 

0900 

6 

1900 

11 

8-21 

1000 

6 

2000 

11 

90S 

1000 

6 

2110 

7.8 

9-06 

1210 

6.4 

2329 

10.3 

O  Night  sample  followed 

8-16,17 

2000 

11 

by  day  sample 

0900 

7 

9-6,7 

2329 

10.3 

1100 

7 

A  Two  night  samples 

8-6,7 

1900 

11 

surrounding  a  day  sample 

0900 

6 

1900 

11 

9-5,6 

2110 

7.8 

1210 

6.4 

2329 

10.3 

TABLE  3-3.  Average  Matt  Balance  of  Ouimette's  Daytime  and  Nighttime 
Fine  Particle  Samples  at  China  Lake. 


Particle 

Average  mass  balance 
ftg/m3  (%  of  accounted  for) 

Day  (7  samples),  % 

Night  (6  samples),  % 

Black  carbon  (soot) 

0.5  (8) 

0.7  (7) 

Organic  C  -  1.5  (total  C  -  soot) 

1.5  (24) 

2.6  (25) 

Sulfate 

2.0  (32) 

2.9  (28) 

Nitrate 

0.1  (2) 

0.2  (2) 

Ammonium 

0.6  (10) 

1.0  (10) 

Crustal  -  ’4(4.2  Si  +  21.6  Fe) 

1.4  (22) 

2.6  (25) 

K  and  Pb 

0.2  (3) 

0.3  (3) 

Unaccounted  for 

5.3 

5.5 
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To  investigate  the  data  for  individual  sampling  days,  Figure  3-11  plots  daytime  values 
(similar  to  8-hour  samples)  versus  daytime  and  nighttime  averages  (similar  to  24-hour  samples)  for 
various  chemical  constituents.  The  plotting  notation  is  explained  in  Table  3-2.  Figure  3-11  reveals 
that  on  individual  days,  the  daytime  data  are  often  much  different  from  the  overall  day  and 
night  averages.  There  are  two  potential  explanations.  First,  the  particle  source  mixtures  may  vary 
greatly  from  day  to  night  on  individual  days  even  though  they  do  not  vary  much  in  an  overall 
average  sense.  This  variation  would  support  the  need  to  consider  shorter  sampling  times  for 
RESOLVE  particulate  measurements.  Alternatively,  the  scatter  in  Figure  3-11  may  represent 
imprecision  in  the  individual  measurements.  Such  imprecision  does  not  appear  unlikely  given  the 
short  sampling  periods  used  by  Ouimette  and  the  subsequent  lack  of  large  deposits  on  the  particle 
filters.  This  would  argue  the  need  for  24-hour  samples  to  obtain  adequate  filter  deposits.  We  are 
not  sure  which  effect  is  most  pronounced.  The  real  answer  may  be  a  combination  of  both  factors, 
although  our  study  of  the  data  leads  us  to  suspect  measurement  imprecision  as  the  cause  of  the 
most  obvious  anomalies. 


3.4  CONCLUSION 

In  summary,  we  have  seen  that  important  diurnal  variations  exist  in  fine  particle  concen¬ 
trations.  However,  the  nephelometry  data  from  the  RESOLVE  program  will  serve  as  a  continuous 
measurement  of  these  variations.  With  respect  to  source  mixtures,  the  high  correlation  between 
8-hour  and  24-hour  bsp  indicates  qualitatively  that  daytime  and  24-hour  particulate  concentrations 
are  not  uncoupled.  Furthermore,  Ouimette’s  data  set,  although  only  for  a  few  days  in  one  season, 
suggests  that  daytime  and  nighttime  particle  compositions  are  the  same  in  an  average  sense. 
Because  of  the  general  implications  of  the  above  analyses,  because  of  procedural  difficulties  and 
extra  expenses  associated  with  short-term  particulate  sampling,  and  because  of  some  interest  in 
nighttime  as  well  as  daytime  visibility  for  DOD  operations,  it  seems  most  prudent  at  this  time  to 
begin  the  routine  RESOLVE  program  with  24-hour  particulate  sampling. 


It  is  recommended,  however,  that  further  measurements  be  taken  and  analyzed  to  check 
this  conclusion.  The  drum  sampler  data  for  the  RESOLVE  program  will  provide  routine  measure¬ 
ments  of  time-resolved  particulate  composition  with  respect  to  soil  elements  and  sulfur.  Further¬ 
more,  as  part  of  the  special  intensive  studies  for  RESOLVE,  the  particulate  samplers  will  be  run 
on  a  daytime  and  nighttime  basis.  By  analyzing  the  drum  sampler  data  and  the  special  study 
data,  it  should  be  possible  to  settle  the  issue  of  8-hour  versus  24-hour  sampling  with  more 
certainty.  If  it  is  found  that,  counter  to  our  current  expectations,  the  particle  source  mixtures 
do  vary  significantly  from  day  to  night,  then  the  sampling  protocol  may  be  changed  later  in  the 
program  to  provide  definitive  data  characterizing  this  effect. 
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POTASSIUM 


4.0  ERROR  ANALYSIS  OF  MASS  BALANCE  AND  PARTICLE 
SCATTERING  BUDGET  MODELS  FOR  RESOLVE 


The  RESOLVE  network  intends  to  quantify  the  contributions  of  particle  emitting  sources 

to  the  extinction  of  light  in  the  atmosphere.  This  quantification  is  to  be  accomplished  via  the 

simultaneous  measurement  of  particulate  chemical  concentrations,  light  scattering,  and  light 

absorption.  The  contribution  of  each  chemical  species  to  the  tight  scattering  is  determined  by  a 
least  squares  solution  to  a  set  of  overdetermined  equations  expressing  the  scattering  coefficient  as 
a  linear  combination  of  the  chemical  concentrations.  This  model  is  termed  the  “particle  scattering 
budget.”  When  the  light  absorption  due  to  elemental  carbon  is  added,  the  model  becomes  an 
“extinction  budget.” 

Source  contributions  to  the  various  chemical  species  are  similarly  determined  by  a 

least  squares  solution  of  simultaneous  equations  expressing  the  chemical  concentrations  as  linear 
sums  of  the  source  contibutions.  This  model  is  termed  the  “mass  balance.” 

The  results  obtained  from  these  calculations  will  have  a  degree  of  uncertainty  arising  from 
two  components:  measurement  uncertainty  and  model  uncertainty.  Measurement  uncertainty  results 
from  the  imperfection  with  which  any  measurement  quantifies  the  true  value  of  a  variable.  It  is 
expressed  as  an  interval  around  a  nominal  value  for  that  variable.  If  these  intervals  are  known  for 
the  measured  values  given  to  a  model,  they  can  be  combined  to  give  an  estimate  of  the  measure¬ 
ment  uncertainty  of  the  model  results. 

Model  uncertainty  result!  from  the  inability  of  the  mathematical  model  to  represent  reality. 
Where  the  discrepancies  between  the  model  and  reality  can  be  identified  and  quantified,  it  is 
possible  to  assign  a  value  to  model  uncertainties  in  various  applications. 

The  treatment  given  here  deals  with  the  random  measurement  uncertainties  of  the  mass 
balance  and  particle  scattering  budget  receptor  models.  Its  objectives  are: 

1.  To  derive  analytical  error  propagation  schemes  for  calculating  the  measurement 
uncertainty  of  results  from  the  mass  balance  and  particle  scattering  budget  models. 

2.  To  calculate  model  results  and  their  measurement  uncertainties  and  verify  them  by 
repeated  calculations  using  randomly  perturbed  input  data. 

3.  To  estimate  the  uncertainties  of  mass  balance  and  scattering  budget  model  predictions 
as  a  function  of  the  uncertainty  level  of  model  input  measurements. 

Owing  to  financial  constraints,  the  scope  of  this  project  was  severely  limited  and  it  cannot 
answer  all  questions  regarding  measurement  and  model  uncertainties  in  RESOLVE.  Several  of  the 
most  important  questions  that  can  be  answered  by  following  this  methodology  will  be  posed  at 
the  end  of  this  treatment. 

4.1  DERIVATION  OF  UNCERTAINTY  FORMULAE 

The  same  solution  method  is  used  for  both  mass  balance  and  particle  scattering  budget 
models. 


NWC  TP  6567 


The  mass  balance  receptor  model  is 


P 

=  2  ^ij  ^j»  *=l>n 
j=l 


(1) 


where 

Cj  =  exact  receptor  concentration  of  species  i 

ay  =  exact  fraction  of  species  i  in  source-type  j 

n  =  total  number  of  species 

p  =  total  number  of  source  types 

The  overbars  are  used  to  indicate  the  value  in  the  absence  of  measurement  uncertainty.  The 
measured  value  will  be  represented  by  the  same  symbol  without  the  overbar.  If  p  is  less  than  n, 
and  Cj  and  ay  are  measured  with  their  corresponding  uncertainties  oq.  and  oa.^  ^  Qf 

equations  (1)  can  be  solved  for  the  most  probable  values  of  Sj  by  minimizing 


(Cj  -  Cj)2 

— - - 

°  C{ 


(ay  -  Ijj)2 

*  .g  ■■  i  - 

o2 


P 


M^i  *  53  »ijSj) 
j*l 


(2) 


where  Xj  represents  a  set  of  LaGrange  multipliers.  The  derivation  is  accomplished  by  the  following 
steps  (References  8  and  9): 

1.  The  derivatives  of  x2  with  respect  to  Cj,  ay,  Xj,  and  Sj  are  set  equal  to  zero. 

2.  Equation  (1)  is  expressed  as  a  Taylor’s  expansion  about  estimates  of  Cj,  ajj,  and  Sj. 

3.  The  expansion  is  combined  with  the  derivatives  acquired  in  step  1  to  express  the  most 
recent  estimate  in  terms  of  the  previous  estimate. 

4.  The  calculated  Sj  values  are  Taylor  expanded  about  estimates  of  the  unknown  true 
value,  S:,  and  the  covariance  matrix  of  the  5j  differences  is  obtained.  The  diagonal 
elements  of  this  matrix  are  the  uncertainties  of  the  source  contributions. 

This  derivation  results  in  the  following  set  of  steps  to  solve  the  mass  balance  equations: 

1.  The  initial  values  of  the  source  contributions  are  set  equal  to  zero. 

*  7  (3) 

is  a  column  vector  containing  the  elements  S^.  The  k  superscript  refers  to  the 
iteration  number. 
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I 


2. 


The  effective  variance  matrix,  V^,  is  calculated 
P 

£  «  &  ♦  £  ^(Sf)2 

j=‘ 


(4) 


VC  is  a  diagonal  matrix  containing  the  t^, ,  the  squared  uncertainties  of  the  ambient 
chemical  concentrations.  is  a  diagonal  matrix  containing  o|..,  the  squared  uncertain¬ 
ties  of  the  source  compositions  for  source  j. 


3.  ~$k+1  =  (AT^)-1  A)-l  aTQ^)'1  T  (5) 

where  A  is  a  nxp  matrix  with  elements  ajj,  the  source  composition  matrix,  and  C  is 

a  column  vector  containing  the  measured  concentrations,  Cj.  Sk+1  is  a  column  vector 
the  (k+l)st  iteration  of  the  source  contributions. 

4.  A  comparison  is  made  between  the  (k-fl)st  and  kth  estimates  of  all  Sj.  If 

SJktl  -  SJK  <0.01,  j=l,p  (6) 

Sf 


then  the  iteration  is  terminated.  If  condition  (6)  is  not  true,  then  steps  2,  3,  and  4 
are  repeated  until  it  is  true. 


5.  The  covariance  matrix  of  the  Sj  is 


Xs  = 


AT  Q*)’1  A 


-1 


(7) 


The  square  roots  of  the  diagonal  elements  of  this  matrix  are  0$., 
of  the  corresponding  Sj.  J 


the  uncertainties 


The  assumptions  of  this  model  and  solution  are  (Reference  8): 

1.  Compositions  of  source  emissions  are  constant  over  the  period  of  ambient  and  source 
sampling. 

2.  Chemical  species  do  not  react  with  each  other;  i.e.,  they  add  linearly. 

3.  All  sources  with  a  potential  for  significantly  contributing  to  the  receptor  have  been 
identified  and  have  had  their  emissions  characterized 

4.  The  number  of  sources  is  less  than  or  equal  to  the  number  of  chemical  species. 

5.  The  source  compositions  are  linearly  independent  of  each  other. 

6.  Measurement  errors  are  random,  uncorrelated,  and  normally  distributed. 
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7.  Measured  source  compositions  closely  approximate  the  true  source  composition,  ay  =  ay. 
Similarly,  the  particle  scattering  budget  model  is 


n-1 

bh  =  oq  +  ^  ai  C|u>  h=l.  in  (8) 

i=l 

where  the  intercept,  Oq,  represents  the  residual  Rayleigh  scattering  in  the  absence  of  particulate 
matter,  and  <*j  is  the  scattering  per  unit  mass  or  scattering  efficiency  of  species  i.  A  nominal 
value  for  Rayleigh  scattering  is  subtracted  from  each  scattering  coefficient  to  obtain  bj,.  Since  this 
value  will  not  be  constant,  a  small  residual  of  non-particle  scattering  is  expected  to  remain,  m  is 
the  number  of  samples. 


If  m  is  greater  than  n,  the  overdetermined  set  of  equation  (8)  can  be  solved  for  the  aj 
using  the  same  procedure  followed  for  the  mass  balance.  The  result  will  be  another  iterative 
solution  in  which 

1.  The  initial  values  of  the  scattering  efficiencies  are  set  equal  to  zero. 

"5°  =  0  (9) 

where  cT  is  a  column  vector  containing  the  elements  j.  The  k  superscript 

refers  to  the  iteration  number. 


n-1 

2-  £  Yq  ($2  (10) 

i-0 

2  2 

where  Vj,  is  a  diagonal  matrix  of  the  and  V(\  is  a  diagonal  matrix  of  the  oq. 

for  each  species  i. 

3.  -5k+1  =  (Cj (V e)~ 1  £)- 1  C^(v£)' *  ?  (11) 

where  C  is  a  mxn  matrix  with  C^q  set  equal  to  unity  to  accommodate  the  constant 
Oq  term.  IT  is  a  mxl  column  vector  with  elements  b^. 

4.  Equations  (10)  and  (11)  are  used  alternately  until 


«r  -  «,k 
~  * 


<0.01,  i=l,  n 


5.  The  oa.  are  the  square  roots  of  the  diagonal  elements  of 


VQ 


CT^-I  £  -I 


(12) 


(13) 
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This  model  is  derived  with  the  following  assumptions  (Reference  10): 

1.  Light  in  a  path  is  scattered  by  only  one  particle  (i.e.,  multiple  scattering  is  ignored). 

2.  All  particles  are  spheres. 

3.  Each  particle  has  a  homogeneous  composition. 

4.  Particle  volume  is  conserved  (i.e„  particles  are  not  aqueous  or  solid  solutions). 

5.  The  index  of  refraction  for  each  particle  is  constant  and  equal  to  the  volume  average 
refractive  index: 

Sr  =  V  £  VjMfj  (14) 

i=l 

where  Mr  is  the  volume  average  refractive  index.  V  is  the  total  particle  volume,  Vj 
is  the  volume  of  particles  of  species  i,  and  Mr,  is  the  refractive  index  of  species  i. 

6.  Mass-normalized  chemical  species  mass  size  distributions  remain  constant  for  all  samples. 

7.  C|y  are  not  significantly  correlated  with  each  other  between  samples. 

8.  Measurement  errors  are  normally  distributed,  random,  and  uncorrelated. 

4.2  MEASUREMENT  UNCERTAINTY  VERIFICATION 

The  analytical  expressions  for  measurement  error  expressed  in  Equations  7  and  13  can  be  veri¬ 
fied  by  simulating  repeated  applications  of  the  measurement  process,  applying  each  model  to  the  set 
of  simulated  measurements,  and  calculating  the  standard  deviation  of  the  model  results.  This 
procedure  is  described  in  detail  by  Watson  (Reference  11).  The  analytical  uncertainties  obtained  from 
Equations  7  and  13  should  not  significantly  differ  from  these  standard  deviations  if  they  are  correct. 

The  advantage  of  a  simulated  test  over  the  use  of  real  measurements  is  that  the  model 
uncertainty  can  be  set  equal  to  zero  and  the  true  measurement  uncertainties  can  be  precisely 
known.  Figures  4-1  and  4-2  illustrate  the  simulation  process  which  is  followed  for  these  tests. 
These  figures  will  be  referred  to  in  the  discussions  below. 

4.2.1  Mass  Balance  Verification 

Source  types  which  could  contribute  to  particle  concentrations  in  the  RESOLVE  monitoring 
region  and  their  compositions  are  presented  in  Table  4-1.  Ammonium  sulfate  is  a  likely  secondary 
product  of  sulfur  dioxide  emissions,  and  much  of  this  gas  is  released  in  the  oil  fields  of  the  lower 
San  Joaquin  Valley.  The  fractional  compositions  of  sulfate,  nitrate,  and  ammonium  in  these  sources 
were  derived  from  their  stoichiometric  ratios.  Windblown  dust  is  considered  to  be  a  major  contributor 
in  a  desert  environment.  Field  burning  is  common  in  California’s  Central  Valley  and  it  is  possible 
that  some  of  its  smoke  could  intrude  into  the  California  desert.  Numerous  roads  passing  through  the 
region  can  add  motor  vehicle  exhaust  to  that  which  is  transported  from  nearby  and  distant  urban 
areas.  The  source  compositions  for  dust,  burning,  and  motor  vehicle  exhaust  in  Table  4-1  were  drawn 
from  Watson  (Reference  11)  and  may  not  be  totally  representative  of  the  California  desert.  These 
source  compositions  are  adequate  for  the  purposes  of  this  study,  however. 
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FIGURE  4-1.  Proem  for  Vorlfyln«  Mm  Balwca  Modal*. 


Case 

Species 

Source  types 

nh4no3 

(NH4)2S04 

Soil 

Burning 

Mot.  veh. 

1 

Mass 

1.000 

1.000 

1.0000 

1.0000 

1.0000 

2 

Organic  C 

0.000 

0.000 

0.0430 

0.4700 

0.5000 

3 

Graphite  C 

0.000 

0.000 

0.0060 

0.0400 

0.0400 

4 

nh4 

0.273 

0.0000 

0.0000 

0.0000 

5 

NO3 

0.000 

0.0000 

0.0200 

0.0090 

6 

0.000 

0.727 

0.0000 

0.0500 

0.0130 

7 

0.000 

0.000 

0.1170 

0.0050 

0.0110 

8 

0.000 

0.000 

0.2540 

0.0050 

00080 

9 

0  000 

0.000 

0.0000 

0.0100 

0.0300 

10 

0.000 

0.000 

0.0100 

0.0650 

0.0070 

11 

Ca 

0.000 

0.000 

0.0090 

0.0090 

0.0120 

12 

Ti 

0.000 

0.000 

0.0080 

0.0007 

0.0000 

13 

V 

0.000 

0.000 

0.0003 

0.0000 

0.0000 

14 

Cr 

0.000 

0.000 

0.0003 

0.0001 

00000 

15 

Vn 

0.000 

0.000 

0.0020 

0.0005 

0.0000 

16 

Fe 

0.000 

0.000 

0.0680 

0.0005 

0.0210 

17 

Ni 

0.000 

0.000 

0.0000 

0.0000 

0.0002 

18 

Cu 

0.000 

0.000 

0.0002 

0.0005 

0.0007 

19 

Zn 

0.000 

0.000 

0.0041 

0.0000 

0.0040 

20 

Br 

0.000 

0.000 

0.0000 

0.0005 

0.0500 

21 

Pb 

0.000 

0.000 

0.0000 

0.0000 

0.2000 

True  source  contributions  of  0.5  pg/W  NH4NO3,  4.5  pg/m^  (NH4)2  SO4.  2  8  pg/m^  soil, 
3.1  pg/m^  burning,  and  0.8  pg/m^  motor  vehicle  exhaust  were  selected  since  these  values  would 
yield  ambient  chemical  concentrations  at  the  high  end  of  the  fine  particle  chemical  concentration 
ranges  found  at  China  Lake  by  Ouimette  and  Flagan  (Reference  10). 

As  Figure  4-1  shows,  these  true  values  of  ajj  and  Sj  were  linearly  combined  to  obtain  true 
ambient  concentrations,  C,.  The  measurement  process  was  simulated  by  adding  a  random  error 
proportional  to  the  uncertainty  of  each  Cj  measurement.  Fifty  sets  of  these  ambient  concen¬ 
trations  were  created  with  uncertainties  of  10,  20,  and  30%  of  the  Cj.  Similarly,  50  corresponding 
sets  of  ajj  measurements  were  produced  for  uncertainties  of  10,  20,  and  30%  of  the  ajj  in  Table 
4-1.  The  mass  balance  model  of  Equations  3  through  7  was  applied  to  each  data  set,  yielding 
50  source  contribution  values  for  each  source  type  and  uncertainty  level.  The  averages,  standard 
deviations  and  ranges  of  the  50  source  contributions  and  their  uncertainties  were  calculated.  These 
results  are  summarized  in  Table  4-2  for  the  three  uncertainty  levels  that  were  examined. 

For  each  source  type,  Table  4-2  compares  values  of  the  following  parameters  for  0,  10,  20, 
and  30%  uncertainty  levels. 

1  50 

Average  Sj  =  <Sj>  =  ^  £  sjh 

h=l 


(15) 
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The  proximity  of  this  avenge  to  the  0%  uncertainty  value  verifies  that  50  trials  provide  a 
valid  statistical  sample  and  that  the  interval  within  which  the  perturbed  values  occur  is  symmetric. 
This  average  is  within  10%  of  the  true  Sj  for  all  cases  except  for  the  30%  uncertainty  level  for 
soil.  Fifty  cases  appear  to  be  adequate  for  statistical  purposes. 

Range  of  Sj.  The  maximum  and  minimum  values  of  the  source  contributions  calculated 
from  each  of  the  50  cases  estimate  the  extremes  that  might  be  encountered  in  a  single  appli¬ 
cation  of  the  model  for  each  uncertainty  level.  These  extremes  differ  substantially  from  the  true 
value,  by  as  much  as  40%,  even  when  the  measurement  uncertainty  is  as  low  as  10%.  These 
extremes  are  within  three  standard  deviations  of  <Sp>  in  every  case,  and  for  30%  uncertainty 
levels  they  are  within  two  standard  deviations. 

I - 50 - 

Standard  deviation  of  Sj  =  2  (Sjh  ■  (16) 

This  standard  deviation  is  one  representation  of  the  uncertainty  of  the  Sj.  If  the  measurements 
could  be  made  50  times,  this  is  the  uncertainty  that  would  be  calculated.  It  can  be  compared 
to  the  uncertainty  calculated  by  the  application  of  the  model  to  a  single  set  of  measurements. 

50 

Average  cSj  =  <oSj>  =  ^  oSjh  07) 

h=l 

The  analytical  uncertainty  was  calculated  by  each  of  the  50  model  runs.  The  average  oSj  is  within 
10%  of  the  standard  deviation  of  Sjk  in  every  case. 

f—50 

Standard  deviation  of  oSj  =  V49  £  (oSjk  -  <oSj»2  (18) 

V  k=l 

When  combined  with  <oSj>,  this  specifies  the  distribution  of  the  oSj  and  shows  the  extent  to 
which  individual  values  of  oSj  will  vary  owing  to  measurement  uncertainty. 

Range  of  oSj.  The  maximum  and  minimum  values  of  the  uncertainties  calculated  from  each 
of  the  trials  show  the  extremes  which  can  be  reached  owing  to  measurement  uncertainties.  The 
maximum  and  minimum  oSj  do  not  always  occur  in  the  same  cases  as  the  maximum  and 
minimum  Sj.  For  10%  uncertainty  levels,  the  extremes  are  within  three  standard  deviations  of 
the  average.  The  extremes  are  often  outside  of  the  three  standard  deviation  interval  for  20% 
and  30%  uncertainty  levels.  For  example,  the  average  soil  uncertainty  at  the  30%  level  is  1.53, 
but  the  upper  extreme  is  four  times  as  large. 

Exact  oSj.  This  is  the  uncertainty  that  is  achieved  when  the  solution  is  performed  with  mass 
balance  assumption  7  met  exactly.  In  most  of  the  cases,  the  errors  introduced  into  the  a^j  values 
cause  deviations  from  this  assumption  with  a  corresponding  over-  or  under-estimate  of  the 
uncertainty. 

Ratio  of  (std.  dev.  Sj)/(exact  oSj).  This  ratio  should  be  close  to  one  if  the  error  propa¬ 
gation  scheme  is  working  perfectly.  This  ratio  also  has  a  statistical  significance  as  a  reduced  chi- 
square  statistic  with  49  degrees  of  freedom  (one  less  than  the  number  of  cases).  With  purely 
random  sampling  from  a  normal  distribution,  as  was  performed  in  this  study,  there  is  a  50% 
chance  that  this  ratio  will  exceed  1.094,  a  20%  chance  that  it  will  exceed  1.163,  and  a  5% 
chance  that  it  will  exceed  1.350.  All  ratios  in  Table  4-2  are  reasonably  obtainable  from  random 
sampling. 
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The  following  conclusions  can  be  drawn  from  the  data  in  Table  4-2: 

1.  The  ability  of  the  mass  balance  model  to  calculate  the  true  source  contributions 
degrades  as  the  uncertainty  of  the  input  data  increases.  The  standard  deviation  of  the  Sj  increases 
nearly  proportionately  with  increases  in  the  measurement  uncertainty.  For  example,  the  standard 
deviation  of  the  50  burning  contributions  increases  from  0.41  to  0.79  to  1.20  as  the  measurement 
uncertainty  increases  from  10  to  20  to  30%. 

2.  The  ability  of  this  model  to  calculate  the  uncertainty  of  the  source  contributions 

degrades  even  more  rapidly  with  increasing  measurement  uncertainty.  For  example,  the  standard 
deviation  of  the  oS;  for  NH4NO3  increases  fourfold  between  the  10  and  20%  cases,  and  another 
fourfold  between  the  20  and  30%  cases.  The  lower  end  of  the  30%  oSj  for  NH4NO3  range 

always  overlaps  the  upper  end  of  the  20%  oSj  range. 

3.  When  measurement  uncertainties  are  high  (i.e.,  20  or  30%)  the  upper  extremes  of  Sj 
are  typically  within  a  factor  of  two  of  the  true  source  contribution.  The  lower  extremes  can 
differ  by  much  more  than  a  factor  of  two.  The  true  source  contribution  is  usually  contained 
within  an  interval  equal  to  twice  the  calculated  uncertainty  for  the  upper  extremes.  It  is  usually 
outside  of  this  interval  for  calculated  uncertainties  associated  with  the  lower  extremes.  It  is 
impossible  to  know  if  one  is  at  these  extremes  in  a  single  application. 

4.  The  majority  of  mass  balance  calculations  yield  uncertainty  estimates  such  that  the 

true  source  contribution  is  within  a  3  oSj  interval,  and  often  within  a  2  oSj  interval,  about  the 

calculated  source  contribution  For  example,  a  burning  contribution  of  1.9  /ig/m3  (one  standard 
deviation  less  than  the  avenge  value  for  the  30%  uncertainty  level)  typically  has  calculated 
uncertainty  of  0.9.  The  true  3.1  pg/m3  value  is  within  1.8  *ig/m3  (2  oSj)  about  1.9  /ig/m3. 

5.  The  exact  oSj  does  not  differ  from  the  standard  deviation  of  oSj  at  a  statistically 

significant  level  yielded  by  a  chi-squared  test.  Other  values  of  oSj,  calculated  using  ay  which 
differ  from  ay  owing  to  measurement  uncertainty,  do  differ  significantly  from  the  standard 
deviation  of  Sj.  The  range  of  oSj  shows  the  extremes  in  every  case. 


4.2.2  Scattering  Budget  Model  Verification 

The  scattering  budget  model  requires  ambient  chemical  concentrations  on  a  number  of 
samples.  Twenty-five  sets  of  concentrations  were  generated  from  the  true  source  contributions 
presented  in  Table  4-1  via  equation  (1)  and  the  source  compositions  in  Table  4-1.  The  values  in 
Table  4-3  were  chosen  to  provide  ambient  chemical  concentrations  at  the  low,  medium,  and  high 
end  of  the  concentrations  measured  at  NWC  by  Ouimette  and  Flagan  (Reference  10). 

These  concentrations  were  combined  into  the  four  categories  of  Ouimette  and  Flagan: 

sulfates  =  SO4  +  NH4 

organics  =  1.2  X  organic  carbon 

crustal  =  14  X  Fe 

unaccounted  -  total  mass  -  sulfates  -  organics  -  crustal 
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TABLE  4-3.  Source  Contribution!  to  Each  Case. 


ii 

12 

13 

14 

L 

15 

16 

17 

18 

19 

20 

a 

21 

22 

23 

(nh4)2so4 

Soil 

0.01 

1.00 

0.90 

0.25 

2.50 

1.80 

0.50 

4.50 

2.80 

0.25 

2.50 

0.90 

0.50 

4.50 

0.90 

0.01 

1.00 

1.80 

0.01 

1.00 

2.80 

0.01 

1.00 

0.90 

0.01 

1.00 

0.90 

0.25 

2.50 

0.90 

0.50 

4.50 

0.90 

0.01 

1.00 

0.90 

0.01 

1.00 

0.90 

0.01 

1.00 

0.90 

0.01 

1.00 

0.90 

0.25 

2.50 

1.80 

0.50 

4.50 

2.80 

0.25 

2.50 

1.80 

0.50 

4.50 

2.80 

0.25 

2.50 

1.80 

050 

4.50 

2.80 

0.50 

1.00 

1.80 

0.01 

4.50 

1.80 

0.50 

1.00 

2.80 

0.01 

4.50 

2.80 

Low  concentrations 

Medium  concentrations 

High  concentrations 

Medium  secondary  particles  dominate 

High  secondary  particles  dominate 

Medium  soil  dominates 

High  soil  dominates 

Medium  burning  dominates 

High  burning  dominates 

Medium  secondary  and  high  burning 

High  secondary  and  high  burning 

Medium  auto  exhaust 

High  auto  exhaust 

Medium  burning  and  auto  exhaust 

High  burning  and  auto  exhaust 

Medium  secondary  and  soil 

High  secondary  and  soil 

Medium  everything  except  burning 

High  everything  except  burning 

Medium  everything  except  auto 

High  everything  except  auto 

High  nitrate,  low  sulfate,  medium  other 

High  sulfate,  low  nitrate,  medium  other 

High  nitrate,  low  sulfate,  high  other 

High  sulfate,  low  nitrate,  high  other 


A  factor  of  14  was  chosen  for  Fe  instead  of  Ouimette  and  Flagan’s  factor  of  21.6  to  make 

the  crustal  components  consistent  with  soil  source  type  of  Table  4-1  (the  inverse  of  0.068  is 

appropriately  14).  The  total  mass  for  each  case  was  set  equal  to  the  sum  of  the  source  contri¬ 
butions  in  each  row  of  Table  4-3. 

Exact  values  for  b^t  were  generated  for  each  of  the  25  sets  of  sulfate,  organics,  crustal, 

and  unaccounted  values  plus  a  small  residual  according  to  the  formula. 

bj,  =  0.01  +  3.20  (sulfates)h  +  2.46  (organics)h 

+  1.42  (crustal)},  +  2.35  (unaccounted)),  (19) 

The  coefficients  of  the  chemical  concentrations  are  the  theoretical  values  calculated  by 
Ouimette  and  Flagan.  A  constant  term  of  0.01  was  chosen  to  be  small  compared  to  the  other 
products.  It  is  of  the  order  of  the  variation  in  the  Rayleigh  scattering  coefficient.  As  Figure  4-2 
shows,  both  the  bj,  and  Cj,j  were  submitted  to  the  same  randomizing  measurement  process  as 
before  with  10,  20,  and  30%  uncertainty  levels.  This  is  analogous  to  having  three  sets  of  50 
particle  samplers  and  three  ncphelometers  at  the  same  location  operating  simultaneously  for  25 
sample  periods.  One  set  of  samplers  and  one  nephelometcr  provide  data  with  10%  precision,  while 
the  other  two  sets  provide  data  with  20  and  30%  precision,  respectively.  The  results  are  presented 
in  Table  4-4  and  have  the  same  meaning  as  those  reported  for  the  mass  balance. 
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TABLE  4-4.  Averages,  Standard  Deviations,  and  Rangat  of  Scattaring  Efficiencia* 
and  Thair  Uncertainties  at  a  Function  of  Uncertainty.  Unaccounted  species 

are  included. 


Species 


1.  Intercept 


2.  Sulfates 


!  3.  Organics 


4.  Crustal 


5.  Unaccounted 


Parameter'1 


Average  Oq 
Range  of  a; 
Std.  dev  aq 
Average  cqjj 
Std.  dev  Cfcj 
Range  of  0^ 
Exact 
Ratio* 


Average  oq 
Range  of  Oq 
Std.  dev  a; 
Average  a q. 
Std.  dev  o^. 
Range  of  0q. 
Exact  OQj 
Ratio 


Average  Oq 
Range  of  Cq 
Std.  riev  oq 
Average  0^ 
Std.  dev  afcj 
Range  of  Ofcq 
Exact 
Ratio 


Average  a; 
Range  of  Qq 
Std.  dev  aq 
Average  Ooq 
Std.  dev  Ofc,. 
Range  of  0^ 
Exact  Cfc^ 
Ratio 


Average  Cq 
Range  of  Cq 
Std.  dev  aq 
Average  o^ 
Std.  dev  Oq. 
Range  of  0^ 
Exact  0^. 
Ratio 


Uncertainty  level 


0% 


0.01 


10% 


0.40 

-2.0  -  0.71 
0.89 
0.84 
0.08 

0.74  -  1.07 
0.89 
1.00 


3.20 


2.46 


1.42 


2.35 


3.03 

2.2  -  3.7 
0.31 
0.33 
0.02 

0.28  -  0.41 
0.35 
0.89 


2.03 

-0.4  -  4.0 
1.07 
1.31 
0.17 

0.95  -  1.55 
•  1.59 

0.67 


1.38 

0.6  -  2.5 
0.42 
0.46 
0.04 

0.38  -  0.55 
0.47 
0.89 


2.95 

0.12  -  6.0 

1.61 

1.96 
0.24 

1.29  -  2.3 
2.38 
0.62 


a  All  units  except  those  of  ratio  are  in  m^/g, 
std.  dev.  Cq 


ratio 


exact  Ofej 
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20% 


0.79 

-3.8  -  4.5 
1.55 
1.49 
0.23 

1.19  -  2.21 
1.78 
0.87 


2.74 

1.1  -  3.9 
0.56 
0.58 
0.06 

0.49  -  0.78 
0.70 
0.80 


1.71 

-1.3  -  5.6 
1.57 
1.84 
0.29 

1.22  -  2.28 
3.19 
0.49 


1.41 

0.0  -  2.5 
0.75 
0.86 
0.13 

0.63  -  1.1 
0.94 
0.80 


3.51 

-1.8  -  8.3 
2.42 

2.76 
0.45 

2.09  -  3.60 

4.76 
0.51 


30% 


1.08 

-4.0  -  5.6 
2.10 
1.96 
0.35 

1.43  -  2.97 
2.67 
0.79 


2.44 

0.40  -  3.8 
0.75 
079 
0.96 

0.65  -  1.97 
1.05 
0.71 


1.68 

-2.1  -  5.1 
1.86 
2.06 
0.34 

1.53  -  2.51 
4.78 
0.39 


1.49 

-0.5  -  3.2 
0.98 
1.17 
0.20 

0.76  -  1.6 
1.40 
0.70 


3.65 
-2.5  -  • 
3.00 
3.13 
0.57 
1.93  -  i 
7.15 
0.42 


Table  44  contains  several  unexpected  results.  First,  the  value  for  the  intercept  is  not  accu¬ 
rately  reproduced  by  the  least  squares  solution.  Its  average  increases  from  a  factor  of  40  to  a 
factor  of  100  over  true  value  as  the  uncertainty  level  increases  from  10  to  30%.  In  each  of  these 
cases,  however,  even  the  minimum  defines  an  interval  about  the  intercept  value  that  contains 
the  true  value,  and  usually  even  surpasses  the  calculated  value  of  the  intercept.  The  interpretation 

can  usually  be  drawn  that  the  intercept  is  not  significantly  different  from  zero  and  should  be 

eliminated  from  the  solution. 

The  average  values  for  each  of  the  <q  do  not  approximate  the  true  values  nearly  as  well  as 
they  did  in  the  case  of  the  mass  balance.  For  example,  the  unaccounted  <a(>  at  the  30% 
uncertainty  level  is  55%  higher  than  its  true  value.  This  means  that  50  samples  are  not  adequate, 
which  is  unlikely  given  the  results  observed  with  the  mass  balance,  or  that  the  variability  about 
the  average  is  not  symmetric.  Many  of  the  low  extreme  values  are  negative.  Though  these  are 

perfectly  reasonable  solutions  to  the  equation,  they  are  physically  impossible. 

Finally,  the  agreement  between  the  exact  and  the  standard  deviation  of  cq  is  not  good. 
Many  of  the  ratios  are  significantly  less  than  1,  some  being  less  than  0.5.  The  uncertainty  calcu¬ 
lation  overpredicts  the  real  variability  due  to  measurement  uncertainty. 

It  is  suspected  that  one  or  more  of  the  model  assumptions  is  violated  by  this  application, 
and  the  most  likely  one  is  assumption  7.  By  definition,  the  unaccounted  variable  contains  a 
linear  sum  of  the  other  variables  and  must  be  correlated  with  them  to  some  extent.  To  test 

this  hypothesis  the  randomized  data  sets  were  regenerated  using  the  equation 

bjj  =  0.01  +  3.20  (sulfates)h  +  2.46  (organics)h  +  1.42  (crustal)j,  (20) 

The  results  are  presented  in  Table  4-5.  There  are  several  similarities  and  differences  with 
respect  to  Table  44.  The  average  intercept  is  still  large  compared  to  the  true  value,  but  it  is 
still  within  the  uncertainty  interval.  The  average  values  of  Oj  more  closely  approximate  the  true 
value  Sj;  the  maximum  differences  are  15%  for  sulfates  at  the  30%  uncertainty  level  in  Table  4-5 
and  is  24%  in  Table  44.  There  are  fewer  negative  extreme  values  for  cq,  and  where  they  exist, 
the  oa.  uncertainty  interval  associated  with  them  includes  values  greater  than  zero.  There  is  close 
agreement  between  the  exact  o <*.  and  the  standard  deviation  of  tq  as  evidenced  by  the  closeness 
of  the  ratios  to  unity.  There  is  probably  a  minor  deviation  from  assumption  7  in  this  case  owing 

to  the  Cju  being  linear  combinations  of  the  same  Sjj,,  but  this  deviation  does  not  appear  to  be 

as  significant  as  it  might  be  in  the  case  of  Table  44. 

Conclusions  similar  to  those  of  the  mass  balance  evaluation  can  be  drawn  from  Tables  44 
and  4-5. 

1.  The  model’s  calculations  of  the  correct  cq  and  oa.  become  less  accurate  as  the  measure¬ 
ment  uncertainty  increases.  In  both  Tables  44  and  45  the  increase  in  the  standard  deviation  of 
cq  is  nearly  linear,  similar  to  the  mass  balance.  For  the  unaccounted  component  in  Table  44, 
the  standard  deviation  of  cq  does  not  increase  as  rapidly  with  increasing  uncertainty;  it  less  than 
doubles  when  the  uncertainty  level  increases  from  10  to  30%. 

2.  The  standard  deviation  of  oa.  also  increases  with  increasing  measurement  uncertainty, 

but  not  so  rapidly  as  it  did  in  the  mass  balance  evaluation.  For  example,  in  Table  44,  oa.  for 

organics  only  doubles  when  the  uncertainty  level  increases  from  10  to  30%.  It  increases  by  a 

factor  of  6  for  the  corresponding  entry  in  Table  4-5. 
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3.  The  upper  extremes  of  S;  calculations  at  30%  uncertainty  levels  always  exceed  the 
true  value  by  a  factor  of  2,  except  for  sulfate,  in  both  Tables  4-4  and  4-5.  The  lower  extremes 
are  not  significantly  different  from  zero,  except  for  sulfate.  This  is  also  the  case  for  uncertainty 
levels  of  20%  when  unaccounted  species  are  included.  There  is  a  reasonable  chance  that  this  model 
may  substantially  over-  or  under-predict  the  true  values  when  measurement  uncertainties  exceed 
20%. 


4.  For  the  majority  of  the  cases  at  all  uncertainty  levels,  the  true  value  falls  within  an 
interval  equal  to  2oa.  or  3oa.  about  the  calculated  value. 

5.  The  exact  uncertainties  calculated  by  the  model  which  includes  unaccounted  species 
overestimate  the  standard  deviation  of  cq  by  a  substantial  amount.  The  exact  do  not  signifi¬ 
cantly  differ  from  the  standard  deviations  of  <*j  when  the  unaccounted  species  do  not  contribute 
to  the  scattering  and  are  not  included  in  the  model. 


4.3  FURTHER  ERROR  ANALYSIS 

This  study  has  dealt  only  with  a  limited  case  and  its  results  may  not  be  applicable  to 
every  case  in  the  RESOLVE  program.  It  has  also  confined  itself  to  measurement  uncertainty  and 
has  not  examined  the  effects  of  model  uncertainty;  i.e.,  the  effects  of  deviations  from  model 
assumptions.  The  methodology  described  here  is  perfectly  applicable  to  the  evaluation  of  model 
as  well  as  measurement  uncertainty. 

As  RESOLVE  progresses,  this  type  of  analysis  should  be  made  an  integral  part  of  the  data 
interpretation  process.  After  the  first  mass  balances  and  extinction  budgets  have  been  performed, 
the  calculated  source  contributions  and  extinction  efficiencies  should  be  used  as  the  true  values 
to  generate  the  randomized  data  sets.  This  will  allow  general  conclusions  to  be  drawn  from  cases 
which  are  more  representative  of  the  actual  situation  than  the  hypothetical  example  used  in  this 
study. 


Known  deviations  from  each  of  the  model  assumptions  should  be  introduced  into  the  simu¬ 
lated  data  generation  scheme.  The  most  important  tests  for  the  mass  balance  model  are 

1.  Test  of  Assumption  5.  The  similarity  of  source  compositions  that  can  still  be  resolved 
by  the  mass  balance  model  can  be  determined  by  making  successive  calculations  in  which  the 
compositions  become  more  equal  to  each  other.  Regional  aerosol  composition  patterns,  which  are 
normally  very  similar,  can  be  evaluated  as  pseudo-sources  to  determine  the  contributions  of  various 
source  regions  to  the  receptors. 

2.  Test  of  Assumption  2.  Nonlinear  interactions,  particularly  those  involved  in  the  forma¬ 
tion  of  secondary  aerosol,  can  be  introduced  into  the  randomized  data  generation  scheme. 
Deviations  from  this  assumption  relate  to  the  accuracy  with  which  primary  and  secondary  sulfate, 
nitrate,  and  organic  material  can  be  apportioned. 

3.  Test  of  Assumption  3.  The  randomized  data  can  be  generated  with  a  greater  or  lesser 
number  of  sources  than  are  included  in  the  solution.  This  will  determine  the  likelihood  of  missing 
a  source  which  is  really  there  or  identifying  contributions  from  a  nonexistent  source. 
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For  the  scattering  budget  model,  the  following  tests  are  important: 

1.  Test  of  Assumption  2.  The  true  scattering  efficiencies  of  randomly  oriented  spherical 
particles  should  be  calculated  and  compared  to  the  spherical  values.  If  the  differences  between 
them  are  smaller  than  the  typical  uncertainties  calculated  by  the  scattering  budget  model,  then 
this  assumption  can  be  relaxed. 

2.  Test  of  Assumption  6.  The  scattering  efficiencies  of  various  size  distributions  should  be 
used  to  generate  the  randomized  data.  The  variability  of  the  model  results  owing  to  this  deviation 
can  then  be  quantified  and  compared  to  the  uncertainties  to  judge  its  significance. 

3.  Test  of  Assumption  7.  Known  correlations  of  increasing  magnitudes  can  be  introduced 
into  the  independent  variables.  This  test  may  resolve  the  curiosities  observed  when  the 
unaccounted  species  were  included  in  the  model. 


4.4  CONCLUSIONS 

Least  squares  solutions  to  the  mass  balance  and  scattering  budget  receptor  models  have  been 
derived  that  provide  uncertainty  estimates  for  model  results.  These  uncertainties  are  functions  of 
the  uncertainties  of  both  the  dependent  and  independent  variables  in  the  models.  Applications 
of  these  models  to  randomized  data  derived  from  known  source  contributions  and  known 
scattering  efficiencies  show  that  the  calculated  uncertainties  are  good  estimates  of  the  real 
uncertainties,  even  when  the  measurement  uncertainty  levels  are  as  high  as  30%.  To  characterize 
and  minimize  the  uncertainty  of  the  interpretive  analysis  results,  the  RESOLVE  program  should 
document  and  minimize  the  measurement  uncertainties  and  maximize  the  number  of  sampling 
cases.  The  latter  would  allow  multiple  applications  of  statistical  models  for  various  data  subjects. 
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Appendix  A 

MONTHLY  DIURNAL  PATTERNS  IN  PARTICLE  SCATTERING 
FOR  FIVE  CALIFORNIA  DESERT  SITES 
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froc.  diff.  from  daily  mean 
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FIGURE  A -3.  (ContdL) 
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FIGURE  A-3,  (Contd.) 
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